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Abstract 

Here we examine the detectability of coUisionless dark matter candi- 
dates that may constitute not all but only a subdominant component of 
galactic cold dark matter. We show that current axion searches are not 
suited for a subdominant component, while direct WIMP searches would 
not be severely affected by the reduced density. In fact, the direct detec- 
tion rates of neutralinos stay almost constant even if neutralinos constitute 
1% of the halo dark matter. Only for lower densities do the rates decrease 
with density. Even neutralinos accounting for only 10~* of the local dark 
halo density are within proposed future discovery limits. We comment 
also on indirect WIMP searches. 

Claims for the need of coUisional cold dark matter [Q as the main form 
of dark matter in the Universe have led us to consider the observability of 
coUisionless cold dark matter (CCDM) when it is merely a subdominant com- 
ponent of the cold dark matter (CDM). Namely, if the previously favored CDM 
candidates, such as axions or Weakly Interacting Massive Particles (WIMPs), 
constitute only a fraction, say 1% or less, of the local dark matter density, would 
these particles still be observable in the current and proposed direct and indirect 
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dark matter searches? This is a vahd question even if non-CCDM is proven not 
to be necessary. In fact there is always the possibiUty of the CDM consisting of 
several populations, the one we are searching for not being the dominant one. 
We could even reverse our question in the following manner. If we see a CDM 
signal in any of our searches, could we be observing a subdominant component 
of the total CDM? 

Naively one may claim that if the local CDM density is 1%, say, of the local 
halo density, the expected rates in CDM detectors, being proportional to the 
local number density, should decrease by the same amount. However, we note 
that a reduction in the relic CDM density implies in general an increase in the 
probability of interaction of CDM with the detector, for example an increase in 
the WIMP-nucleus cross section or an increase in the axion-photon coupling 
constant. Since the detection rate depends on the product of the interaction 
probability and the local CDM density, the increase in interaction probability 
may compensate the decrease in CDM density, and the detection rate would 
remain unchanged. 

For axions, this argument is new; for WIMPs, it is not. It has been men- 
tioned implicitly or explicitly in many papers on WIMP detectability since the 
inception of the subject[^. It is timely, we believe, to pinpoint, emphasize and 
update this argument, because it clearly points to the value of continuing WIMP 
searches even if WIMPs constitute only a small fraction of the dark matter. 

We now present arguments that the compensation between interaction prob- 
ability and local density occurs for axions and WIMPs, and point out some 
exceptions. 

Unless there is segregation for different types of dark matter, the ratio of 
CCDM to total DM should be the same locally in the Galaxy and globally in 
the whole Universe. Thus in the following we assume that the local fraction of 
CCDM fee DM is related to the CCDM relic density ^Iccdm through 

r PCCDM ^CCDM 

J CCDM — — , (i-) 

Plocal i 'DAf 

where pccDM is the local density of a particular CCDM candidate, piocai — 0.3 
GeV/cm^ is the local halo density (at the location of the Earth), ^Iccdm is the 
relic density of our particular CCDM candidate, and i^oM — 0.3 is the total 
contribution of DM to the total energy density of the Universe. 

Because the relic density of axions is directly related to its mass, and axion 
searches are tuned to the axion mass, current searches are not suited to look 
for a subdominant axion component. The axion relic density is directly related 
to its mass ma- The usual relation (which has its caveats, see for example 
and references therein) between the axion relic density and its mass rua is, for 
a QCD constant of 200 MeV, 

0.6 X lO^^eV 

where h is the reduced Hubble constant, h ~ 0.7. A dominant component of 
axions with fla = 0.3 corresponds, according to this relation, to ma = 3 x 
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10~^eV. Thus, we could decrease the density at most to V,a = 0.003, so that 
axions contribute 1% of the total DM density, before encountering the upper 
bound of 3 X 10^'^eV on the axion mass derived from the observed duration 
of the Supernova 1987A neutrino signal (and other bounds which exclude all 
heavier axions, see for example and references therein). 

The power P from axion to photon conversion in an electromagnetic cavity 
used for axion dark matter searches is proportional to the product PaTUa of the 
local axion density and the axion mass d). In absence of segregation, eq. (0) 
shows that the power is also proportional to flama, which using eq. @) for the 
axion relic density gives 

pcx^]y^ (3) 

that is the power is proportional only to the l/7th power of the axion relic 
density. For a decrease in fla by a factor of 100, the power decreases only by 
a factor of 2. Of course, because the axion mass has shifted to keep relation 

valid, this power is now at a frequency which is 500 times larger and one 
would need resonant cavities consequently smaller. The limiting factor of axion 
dark matter searches with electromagnetic cavities is not the axion to photon 
conversion power, but the size of the necessary cavities. 

The relic density of WIMPs fi^ is determined by their annihilation cross 
section aa by the relation 



n^h' ^ ^ , (4) 



1 X IQ- 

where (aav) is the thermal average of the annihilation cross section times the 
relative velocity of the WIMPs at freeze-out. A reduction in the relic WIMP 
density requires an increase in their annihilation cross section in the early Uni- 
verse. This increase is often associated with an increase in the scattering cross 
section as of WIMPs off atomic nuclei. Since the interaction rate in detec- 
tors depends on the product OsPxi if the scattering cross section increases as 
much as the annihilation cross section, the rate would be unchanged even if 

has decreased. Concerning indirect detection, the flux of rare cosmic rays 
and of gamma-rays produced in halo annihilations depends on the product of 
the square of the density and the annihilation cross section into a particular 
channel, OaP^ ■ Thus, even if an increase in the cross section would compensate 
the decrease in one of the powers of the density, the fluxes would still decrease 
linearly with the halo WIMP density. However, the intensity of the high-energy 
neutrino emission from the Sun and the Earth would in many cases decrease 
only slightly, because, to the extent that capture and annihilation of WIMPs 
in the Sun and the Earth have the time to equilibrate, the neutrino intensity 
depends only on the capture rate which in turn depends on the product (JsPx- 

We can understand the relation between the scattering and annihilation cross 
sections (jg and (Tq as follows. The scattering cross section of a WIMP of mass 

with a nucleus of mass m^r is of the form 

" :\As\' , (5) 



(m^ -I- m,N 



|2 I 
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where Ag is a reduced amplitude which depends on the dynamics of the coUision. 
The annihilation cross section of WIMPs into light particles is 



<7„~ Af„m^|A„|^ (6) 

where Aa is the corresponding reduced amplitude and Na is the number of 
annihilation channels. In the case of interactions of weak order, the amplitudes 
are of the order, 

\As\^^A^^^, (7) 

where a is a coupling constant of weak order a ~ 10~^, M is a mass of the 
particles mediating the interaction, typically M ~ lOOGeV and A is the atomic 
number of the interacting nucleus. Our expression for the scattering amplitude 
includes the nuclear coherent enhancement factor valid for spin-independent 
scattering; for spin-dependent scattering the factor A^ should be dropped. Also, 
our expression for the annihilation cross section is valid for m^^ < M, while in 
the opposite range, > M, we expect aa — Nam~^. 

The simplest case to consider is that of WIMPs lighter than the nuclei they 
interact with. Prom the above equations it is obvious that for these WIMPs 

the ratio of cross sections is approximately constant. In fact, provided the main 
annihilation channel is into fermions, quarks in particular, crossing arguments 
insure that the reduced amplitudes of annihilation and scattering with nucleons 
are similar. 

Heavier WIMPs may have other annihilation channels, such as Higgs bosons 
or vector boson pairs. The crossing argument then does not apply and we don't 
expect the scattering amplitude to grow as much as the annihilation amplitude. 
Moreover, for WIMPs heavier than the nuclei they scatter from, the scatter- 
ing cross section becomes largely independent of the WIMP mass, while the 
annihilation cross section always depends on m,^. In this case, while the annihi- 
lation cross section could be made larger by considering lighter (if > M) or 
heavier (if < M) WIMPs, the scattering cross section would remain largely 
unchanged. 

Therefore, for relatively light WIMPs, and to a lesser extent for heavy 
WIMPs, we expect the scattering cross section to grow by the same factor 
VIdm/^x annihilation cross section needs to grow to reduce the local CDM 
density by n^/floM- So the rate, which is proportional to the product of the 
local CDM density and the scattering cross section, remains unchanged. 

This argument ceases to be applicable at some small enough WIMP densi- 
ties, because the necessary increase in cross sections is due to larger couplings 
and/or smaller mediator masses, which, at some point, encounter accelerator 
limits which exclude the model. In fact Fig. 1 (described below) shows that for 
neutralinos constituting 10% of the halo or more the direct detection rates are 
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largely maintained (as evidenced by the behavior of the envelope of the highest 
rates), and for densities as low as 1% of the halo density, the highest rates only 
decrease by a factor of about three, showing that there is compensation in the 
interaction rates while densities decrease by a factor of up to 100. As mentioned, 
the compensation ceases to work for smaller densities, and for these (as can be 
seen in Figs. 1 and 2) the envelope of highest rates decreases linearly with the 
density. 

To substantiate the general arguments presented so far, we have analyzed the 
concrete case of the lightest neutralino in usual variations of the Minimal Super- 
symmetric Standard model. We used a table of models allowed by all accelerator 
limits, produced with the DarkSUSY code Q over the last few years for other 
purposes, i.e. having in mind other issues which were addressed in the papers of 
Ref [H for which the models were originally computed. We have, therefore, not 
done any particular sampling of the models to favor lower densities and higher 
detection rates. We restricted our attention to models with fi;^ < ^Idm = 0.3 
(r2ji,/i^ < 0.15) for which we found about 45,000 points in parameter space. 
For these models, using the spin-dependent and spin-independent neutralino- 
nucleon cross sections provided in the table, we computed the integrated inter- 
action rates on Ge, following L. Bergstrom and P. Gondolo in ref. Q. We plot 
the resulting integrated rates (in units of events per kg-day) in the first two 
figures of this paper. 

Figs. 1 and 2 show the expected integrated rates in Ge detectors as function 
of the lightest neutralino relic density. Fig. 1 shows only a part of Fig. 2 (the 
part with the highest rates and densities) displaying the original points in the 
table of models. Fig. 2 shows the whole range of densities (which reach up to 
fl^h^ ~ 10^®) using a regular grid of points covering the region with models. 

In Fig. 1 the change of the slope of the envelope of the points with maximal 
rate as the density diminishes is clearly evident. There is approximately no 
change in maximal rates in the first decade of decrease of density, from fl^h^ = 
0.15 (for which neutrahnos constitute the whole halo, fccDM = 1) to 0.015 (for 
which neutralinos constitute 10% of the halo, fccDM = 0-1). There is only 
about a factor of 3 decrease in the next decade, from U^h"^ — 0.015 to 0.0015 
ifcCDM from 0.1 to 0.01). For smaller densities the slope of the envelope clearly 
changes, and as evidenced by Fig. 3, the maximal rates decrease linearly with 
up to the smallest densities. Some of the points shown in Figs. 1 and 
2, mostly among with the smallest densities in Fig. 3, should correspond to 
resonances in the annihilation cross section. 

The compensation in the rates can be largely understood just by looking 
at the spin-independent neutralino-proton cross section c^-p as a function of 
the lightest neutralino relic density, shown in Fig. 3, again with a regular grid 
showing the allowed region where points were found. Also from this figure, 
looking at the envelope of the highest cross sections, it is evident that for Jl^/i^ 
decreasing from 0.15 to 0.0015, i.e. in the first two decades of decrease in 
neutralino density, a^-p increases with decreasing densities; this leads to a 
compensation in the direct rates. On the other hand, for smaller densities, <Tx-p 
is about constant or decreases slightly with decreasing densities; this effect is 
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due to accelerator bounds. 

Since experimental upper bounds and discovery regions are at present given 
in terms of cr^-p, Fig. 3 shows the approximate level of the claimed signal and 
present bounds (by the DAMA, CDMS, COSME-IGEX, and Heidelberg-Moscow 
collaborations [^) and conceivable future discovery level (by the GENIUS pro- 
posal 1^]) which are of order 10~^ pbarns and 10~^ pbarns, respectively, for 
neutralinos which account for the whole local halo density, i.e. with fcCDM=^- 
(These values depend on the neutralino mass, but to simplify the presentation 
we only take the most conservative bounds in our range of masses.) In our 
case these values must be understood as levels of /ccdmCTx-pj which are shown 
in Fig. 3 (with short-dashed and long dashed lines respectively). The present 
level of discovery lightly touches the boundary of the highest rates for densities 
reduced by up to a factor of about 10. This suggests the possibility that the 
DAMA claimed signal may correspond to subdominant neutralinos. It is very 
interesting to see that many models of subdominant neutralinos even with 10~^ 
of the total dark matter density, enter in the discovery limit proposed by Genius. 

In conclusion, the main point of this paper is that the direct detection rates 
of neutralinos remain about constant for neutralino densities between 100% and 
1% of the halo dark matter and only decrease linearly with the density for lower 
densities. Thus if a signal is found in direct detection experiments the question 
of which component of dark matter was found, the primary or a sub-dominant 
one, may remain open. We also note that neutralinos with density as small as 
10~^ of the local dark halo density are within the discovery limits of proposed 
future experiments. 
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Figure Captions 



Fig. 1 Integrated interaction rates of neutralinos in Ge detectors (computed 
as in L. Bergstrom and P. Gondolo Ref. [^) in units of events per kg-day, 
as function of the ncutralino relic density, for il.y,h^ < 0.15. Each point 
represents an actual model. 

Fig. 2 Integrated interaction rates of neutralinos on Ge extended to the whole 
range of densities. A regular grid of points shows the region covered with 
models. 

Fig. 3 Spin-independent neutralino-proton cross section function of 

the lightest ncutralino relic density. As in Fig. 2, a regular grid of points 
shows the region where models were found. The short-dashed and long 
dashed lines of fccDMO'x-p~^^~^ pbarns and 10~^ pbarns show the ap- 
proximate level of DAMA claimed signal and the current bounds, and the 
conceivable future discovery level, respectively {fee dm is the fraction of 
the local halo density consisting of neutralinos). 
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Figure 1. 



9 



Duda, Gelmini, Gondolo - 2001 



102 
101 j- 

^ 10-1 

^ 10-2 ^ 
> 10-3 
§ 10-4 L_ . 

J, 10-5 J 

^ 10-^ I 

£ 10-^ ^ 
S' 10-8 L_ 

10-9 L 



10 



Figure 2. 
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Figure 3. 
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